For PHS 7045
The University of Utah
2024-12-12
Hill-of-vision: a graphical representation of retinal sensitivity, “peaks” represent areas of high sensitivity.
We would like to use this hill-of-vision map to calculate standard sensitivity loss
Conservative method: spatial interpolation
Spatial interpolation estimates sensitivity for a single virtual patient’s testing results using Kriging.
However, this approach may lose some population-level information and is harder to interpret.
An interpretable age-adjusted prediction model would be much better!
Works with multiple exam types (e.g., Mesopic, Cyan, Red, and Cyan-Red difference).
Provides predicted sensitivity values and visualizations with parallel computing.
Kriging uses the best linear unbiased prediction based on both the spatial structure and the measurements, so it’s not guaranteed to be the best prediction, especially in areas with sparse test points.
| ExamType | Model | RMSE | PartialSill | CorrDistance |
|---|---|---|---|---|
| Mesopic | Spherical | 0.3715 | 19 | 6 |
| Dark-adapted Cyan | Spherical | 0.9294 | 24 | 6 |
| Dark-adapted Red | Spherical | 0.5114 | 22 | 10 |
| Cyan-Red difference | Stable | 0.9096 | 22 | 8 |
Objective: Predict retinal sensitivity based on patient-specific characteristics like age and test location.
Models Used:
Linear Mixed Model (LMM), with age and eccentricity.
Bayesian Quantile Regression (BQR), including an isotropic smooth (thin plate regression spline) for the x and y coordinate, a spline for the age and a tensor product smooth for age and eccentricity.
Random Forest (RF), using the variables: participant age, eye, eccentricity and angular distance from the fovea, and spatial coordinates (x,y).
Performance comparison: mean absolute error (MAE) and mean absolute calibration error (MACE), with patient-wise cross-validation and conformal prediction frameworks.
SensForFit()PredictCompare()
PredictNormal_lmm(), PredictNormal_bqr(), PredictNormal_rf()PredictSensitivity()VisualRetina()data("refMes")
## data preparation
refMes <- SensForFit(dt = refMes, examcol = "Examtype", idcol = "Patient", agecol = "Age", senscol = "MeanSens", k = 10)
## model comparison
(res.tab <- PredictCompare(dt = refMes, exam = "Mesopic", CalibSplit = 0.2, coverage = 0.95))
## prediction for a future patient aged 77 (e.g. bqr is the best)
pred_BQR_77 <- PredictSensitivity(model = "BQR", age = 77, dt = refMes)
## hill-of-vision for new patient
VisualRetina(pred_sens = pred_BQR_77, exam = "Mesopic")| Model | MAE | MACE | |
|---|---|---|---|
| lmm_results | LMM | 1.43 | 0.03 |
| bqr_results | BQR | 1.70 | 0.02 |
| rf_results | RF | 1.53 | 0.03 |
Other Models vs. Interpolation
Models provide insights into age and test-location effects on sensitivity.
No need to predict complete hill-of-vision.
Models can account for variability across patients and test sites.
Further exploration is needed to develop a robust conformal prediction framework, especially for machine learning methods.
[1] Pfau, M., Jolly, J. K., Charng, J., von der Emde, L., Müller, P. L., Ansari, G., Pfau, K., Chen, F. K., & Wu, Z. (2024). Multicenter normative data for mesopic microperimetry. Investigative Ophthalmology & Visual Science, 65(12), Article 27. https://doi.org/10.1167/iovs.65.12.27
[2] Pfau, M., Müller, P. L., von der Emde, L., Lindner, M., Möller, P. T., Fleckenstein, M., Holz, F. G., & Schmitz-Valckenberg, S. (2020). Mesopic and dark-adapted two-color fundus-controlled perimetry in geographic atrophy secondary to age-related macular degeneration. Retina, 40(1), 169–180. https://doi.org/10.1097/IAE.0000000000002337
[3] Angelopoulos, A. N., & Bates, S. (2021). A Gentle Introduction to Conformal Prediction and Distribution-Free Uncertainty Quantification. arXiv. https://doi.org/10.48550/arXiv.2107.07511
[4] Johansson, U., Boström, H., Löfström, T. et al. Regression conformal prediction with random forests. Mach Learn 97, 155–176 (2014). https://doi.org/10.1007/s10994-014-5453-0
Yingjia Wei